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We derive the mean-field theory of disordered block-copolymers composed of two monomeric 
species, combining Edwards' functional method with the replica technique of disordered systems. 
In the absence of disorder we recover the canonical self-consistent field theory of inhomogeneous 
polymers. In the presence of sequence disorder the theory can be regarded as a comprehensive, novel 
self-consistent treatment of copolymer melts, unifying the weak- and strong-inhomogeneous regimes. 
In particular, we study the stability of the microphase separation transition in a melt of diblock- 
copolymers, composed of A'^ monomers and equal A- B-volume fraction, against the disruptive effect 
of disorder. We obtain a phase diagram in terms of the relevant parameters, namely the rescaled 
Flory-Huggins parameter a-nd disorder strength p. 
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The interplay between topological and energetic effects 
in linear fieteropolymeric systems is a central problem in 
soft-condensed matter and biophysics. Heteropolymers 
composed of two monomeric species include for example 
block-copolymers, alternating {AB)n multiblock copoly- 
mers and random copolymers [lIQQ. Polymer chains 
composed of several monomeric species include random 
heteropolymers, biopolymers as well as proteins, that are 
comprised of twenty amino-acids [3| [s] ■ 

Block-copolymer melts undergo an order-disorder 
transition at which spatially periodic, coexisting mi- 
crodomains of similar chemical composition form. In 
particular, mean-field theory [7] [s] predicts that a melt 
of diblock-copolymers with an equal A- and B-volume 
fraction undergoes a second order phase transition at 
which a lamellar phase (L) occurs. For different values 
of the A- and _B-volume fraction, diblock-copolymers are 
known to present in addition other morphologies, char- 
acterized by periodic two- and three-dimensional spacial 
orderings, e.g. body-centered-cubic (bcc), hexagonally 
ordered cylindrical (hex) and bicontinous gyroid (G) [9i] 
and mean-field predicts all transitions to be first order 
in this case [lOj . Different microphase orderings, induced 
by the appropriate choice of molecular weight and chain 
architecture during synthesis, have different thermal and 
mechanical properties. Consequently, block-copolymers 
find several applications in surfactant systems, as well as 
microelectronic systems and nanotechnology. 

The phenomenon of microphase separation in systems 
of block-copolymer melts have been explained, within the 
weak-segregation approximation, more than thirty years 
ago [1]. A self-consistent, mean- field theory of block- 
copolymer systems for the gaussian chain model, be- 
yond the weak-segregation approximation, has been for- 
mulated [H and, despite the high computational cost 
required to solve the corresponding equations, it repre- 
sents to date the leading method to investigate quanti- 
tatively block-copolymer morphologies. The phase di- 



agram in terms of the Flory-Huggins parameter x^ the 
degree of polymerization N and A-block volume fraction 
/ has been obtained and represents one of the crowning 
achievements of self-consistent field theory [l^j . 
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FIG. 1. The phase diagram of a disordered diblock-copolymer 
melt with equal A- _B- volume fraction / = 1/2, and disorder 
strength p, in terms of the relevant parameters, namely the 
rescaled Flory-Huggins parameter and disorder strength, 
p. The solid line represents the second-order phase boundary 
between the lamellar (L) and disordered (D) phases. 

In this letter we study the effect of quenched disorder 
on the equilibrium properties of block-copolymers 
The equilibrium mean-field theory here introduced does 
not describe the effects of self -generated disorder, that 
plays an important role in dynamical mean-field theories 
of polymer melts and glassy systems [l^l [l3 • 

The study of condensed-matter systems that possess 
quenched disorder represents a relatively old, yet 
rapidly evolving frontier of modern science. The mean- 
field theory of spin-glasses 15| 16[ has been extended to a 
variety of different problems, including the one of chem- 
ically disordered polymeric systems, firstly at the phe- 
nomenological level 17|, and later from the microscopic 
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point of view [18| . Quenched disorder in polymeric sys- 
tems may occur to be fixed i n sp ace, when studying the 
behavior of macromolecules [ij, manifolds [2^ and di- 
rected polymers in random media 2lj or might instead 
characterize the sequence of monomeric units. 

The equilibrium properties of random AB copolymer 
melts have been studied in the weak-segregation approx- 
imation, by means of the replica [22| and cumulant ex- 
pansion methods [2^. In these two studies, the effect 
of correlated disorder and the presence of finite size A,B 
multiblocks have been considered in detail. An equivalent 
expression for the Landau-type density functional free- 
energy has been obtained with both methods and the cor- 
responding phase diagram, in the weak-segregation limit, 
has been discussed. Related results have been obtained 
earlier, within the ground-state dominance approxima- 
tion [18[, in the case of uncorrelated sequence disorder. 



We derive the self-consistent theory of disordered het- 
eropolymers, that includes both the weak- and strong- 
segregation regimes (2^ 24 1 , and apply it to the particu- 



lar problem of disordered AB diblock-copolymers. 

We consider the canonical formulation pj for a melt 
of Up block-copolymer monodisperse chains composed of 
N segments, where the index i = 1, • • • , rip labels each 
chain, where chains are characterized by A-block volume 
fraction /, where the total density of the system is po, 
where {ri(s)} are the monomer cartesian positions and 
where 0i{s) = ±1 are the occupation variables for seg- 
ments A and B respectively, that occur to be disrupted 
with probability pa,b, in each block independently. The 
total density for segments of both type A and B and the 
relative density of the two species are defined as 

N 

•^W = ^E dsS{r-r,{s)) 

N 

^W = ^E / ds9.{s)S{r~r,{s)), (1) 

and depend on the quenched random occupation vari- 
ables {9i}, where each specific realization encode the 
presence of uncorrelated sequence disorder within each 
block The model here introduced relates, without 

loss of continuity, the block-copolymer problem in the 
absence of disorder {pa,b = 0), with the problem of ran- 
dom AB copolymers {pa,b = 1/2)- The hamiltonian of 
the system is given by 

+ XPo f dr[m^{r)-4>^{r)], (2) 



and where aN^^^ is the statistical segment length of the 
Gaussian chain model. 

We considered the case of uncorrelated disorder, so 
that the random occupation variables obey a binomial 
distribution in each block independently, with variance 
{di{s)6j{s')) = 4:6ijS{s - s')p{l - p), where p = pa,b for 
s < f and s > f respectively, and where we assume 
that monomers A and B are represented by statistical 
segments with the same size aN^^^. The order param- 
eter M is defined as the amplitude of the relative den- 
sity of the two monomeric species, m(r). According to 
the Markov model of random copolymerization 23| , each 
block is characterized by its length and disorder strength, 
chemically wired during synthesis, i.e. an average copoly- 
mer composition pa,b (or mole fraction of monomers of 
type A in block B and viceversa). 
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where x is the Flory-Huggins parameter, measuring the 
relative strength of the interaction between similar and 
dissimilar monomers, where we work in units of fc^T, 
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FIG. 2. The order parameter M obtained from self-consistent 
field theory, for different values of the rescaled Flory-Huggins 
parameter and disorder strength. In the first plot (top) 
we show the order parameter M as a function of disorder 
strength, for increasing values of the rescaled Flory-Huggins 
parameter = 20, 25, 40, 80, 100 respectively. In the sec- 
ond plot (bottom) we show the order parameter Ai as a 
function of x^^ for increasing values of disorder strength 
p = 0,0.05,0.1,0.15, 0.20. 

The two microscopic density operators in equation ([T]) 
can be eliminated from the partition function, introduc- 
ing the four, conjugate fields (j){r), m{r), w{r) and v{r). 
The incompressibility condition of the melt is enforced 
introducing the total chemical potential ^(r). We aver- 
age the partition function over the uncorrelated distri- 
bution of quenched random occupation variables in each 
block independently {Oi{s)) = ±{1 — 2pa,b), according to 
the self-averaging hypothesis. Differently from the prob- 
lem of self-interacting heteropolymers composed of many 
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monomeric species j5| , where a replica symmetry broken 
scheme has to be introduced @, the eigenvalue analysis 
of the AB copolymer free-energy shows that the replica 
symmetric solution is stable [18j, so that the analytic 
continuation to the n = limit is a simple task [23j . 
The self-consistent saddle point equations are given by 
w{r) = 2xN(I>{y) + ^(r) and v{y) = ~2xNm{r). The 
conjugated density fields are given by 



rn(r) 



1 

2Q 
1 

2Q 



dse(s)[l - 2X(r, s)]q{v, s)q\v, s) 



dsq{v,s)q\v, s), 



(3) 



where the single-chain occupation variable, having per- 
formed the average over disorder, is simply given by 
8(s) = +1 for s < f and 6(s) = -1 for s > f and 
where the single-chain polymer partition function is de- 
fined as Q = / drq{r, s)q^ (r, s). The two conjugated den- 
sity fields in equation ([3]) are related to the forward prop- 
agator q{r, s) via the modified diffusion equation, 



(4) 



and a similar modified diffusion equation applies for the 
backward propagator gt(r, s) . The self-consistent mean- 
field, given by >V(r,s) = iw(r) + ^v{r)[l - 2P(r,s)], 
has two contributions, associated to the two conjugated 
density fields w{r) and t;(r) above and depends on the 
disorder strength pa,b via 



2?(r,s) = In [l +p(s)(: 



l)]/e(r,s) 



X(r,s) =p(s)e^('-'^V[l+P(s)(e^''''^ - l)]- (5) 

The two functions in equation ([5]) are obtained averag- 
ing with respect to the binomial form of the disorder 
each term in the partition function and rcsumming the 
associated series and depend on the average, auxiliary 
mean-field C/(r, s), containing four-propagator contribu- 
tions 0), defined by 

g{r,s)=e{s)[v{r)-^l dr'v{r')q{r',s)qHr\s)]. (6) 

In the absence of disorder, the saddle-point equations 
Q-dSl) reduce to canonical self-consistent field theory 
|10| . The modified diffusion equation (U) has been de- 
rived, within a replica symmetric ansatz for the many- 
body propagator q{ri, ■■■ , r„, s) = Yia^i 9(ra, «), follow- 
ing the principles of time-dependent mean-field Hartree- 
approximation ^2^. 

We consider a system of block-copolymers with equal 
A-S-volume fraction and the particular case of a sym- 
metric disorder strength within the two blocks, pa.b — P- 
To solve numerically the associated saddle-point equa- 
tions ©-dH) we used the real-space, one-dimensional 
Crank-Nicholson algorithm, in combination with a sim- 
ple mixing technique. For any given value of the rescaled 
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FIG. 3. Domain spacing d as a function of the rescaled Flory- 
Huggins parameter for different values of the disorder 

strength, p = 0, 0.05, 0.1, 0.15, 0.20. The dashed line indicates 
the asymptotic scaling form predicted by strong-segregation 
theory, d tx (xN)^^^ UM- the order-disorder transition all 
curves, corresponding to different values of disorder strength, 
terminate at the same value d* ~ 1.32, as expected along 
the second order phase boundary separating the lamellar (L) 
and disordered phase (D) and indicated by the filled, black 
squares. 



Flory-Huggins parameter we find the optimal value of 
the domain spacing d that minimizes, according to a 
standard one-dimensional minimization routine, the free- 
energy for different given values of disorder strength and 
rescaled Flory-Huggins parameter x^- Typically we 
used 10'^ points along both the spacial and temporal 
grids, though higher values 5 x 10'^ have been used to 
check accuracy of our results in the strong-segregation 
regime. 

In Fig. [T]we present the phase diagram of disordered 
block-copolymers with A-block volume fraction / = 1/2. 
The solid thick line represents the second order phase 
boundary between the lamellar (L) and disordered (D) 
phases. As expected, lamellar order (L) is inhibited by 
the presence of disorder, but persists in the low temper- 
ature region, at large values of the disorder strength. 

In Fig. [2]we show the values of the order parameter M, 
for different values of disorder strength (top panel) and 
rescaled Flory-Huggins parameter xTV (bottom panel), 
corresponding to horizontal and vertical scans in the 
phase diagram of Fig. [TJ The order parameter M is 
observed to depend on disorder. The location of the 
second-order phase boundary between the lamellar (L) 
and disordered (D) phases in the phase diagram Fig. [T] 
corresponds to the onset of points where the order pa- 
rameter vanishes, as in Fig. [5] The expected asymp- 
totic form M = (1 — 2p), shown as a dashed line in the 
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top panel, is found to agree very well with the numer- 
ical results obtained within self-consistent field theory, 
in the intermediate segregation region ^ 60 and at 
small values of disorder strength p < 0.2, while devia- 
tions form the strong-stretching results are described by 
self-consistent field theory, as in Fig. O The asymptotic 
form we derived above for the order parameter M, scales 
as the effective block-size in the presence of uncorrelated 
disorder and numerical self-consistent field theory is ob- 
served to approach exactly this simple form. 

In Fig. [3] we plot the domain spacing d, character- 
izing lamellar order (L), as a function of the rescaled 
Flory-Huggins parameter x^- Domain spacing is ob- 
served to depend on the disorder strength, but not at 
the order-disorder transition, i.e. along the second-order 
phase boundary in the phase diagram, and as shown by 
the black filled squared boxes. In the strong-segregation 
limit the domain spacing of pure, non-disordered block- 
copolymers approaches the asymptotic form d cx {xNY^^ 
[24| . In fact, we observe that the domain spacing of dis- 
ordered block-copolymer systems appears to cross-over 
to the same universal asymptotic form, shown by the 
dashed line. The calculation of a global phase diagram, 
in the space of disorder strength p, block A volume frac- 
tion / and rescaled Flory-Huggins parameter xiV can be 
addressed within the present theory, and the issue of the 
stability of other block-copolymer morphologies, beyond 
the one-dimensional lamellar (L) ordering here discussed, 
and the nature of the phase transition within the mean- 
field approximation and in the presence of quenched dis- 
order is an open interesting problem. 

The important role of fiuctuation effects, in the frame- 
work of one- loop Hartree approximation [i^ , as well 
as momentum space renormalization-group theory [27| . 
have been considered in the absence of disorder and cor- 
rections to the mean-field phase diagram of pure block- 
copolymer melts have been discussed [1^ [20|. The ef- 
fect of thermal fiuctuations in random copolymer sys- 
tems have been discussed, within the weak-segregation 
approximation [l8| . A model of block-copolymer systems 
with equal A- B- volume fraction / = 1/2, in the presence 
of gaussian distributed Flory-Huggins interaction values, 
characterized by an average x disorder strength Ax, 
and the corresponding phase diagram have been dis- 
cussed recently within renormalization-group theory [s^] ■ 
It has been argued that the presence of disorder changes 
the nature of the transition at a critical value of disorder 
strength and lamellar order becomes unstable at values 
of Ax/Xs — 0.1, where Xs — 10.495 is the order-disorder 
transition in the absence of sequence disorder. 

We observed that, within the mean-field theory ap- 
proximation, lamellar ordering persists to large values of 
disorder strength, for any value of the Flory-Huggins pa- 
rameter our algorithm can reach, and we in fact expect 
lamellar ordering to persist at low enough temperatures, 
for any value of disorder strength. 
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